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d Abstract 

A new classification scheme for pairs of real numbers is given, generalizing work 
from which in turn was motivated by ideas from statistical mechanics in general 
I and work of Knauf [17] and Fiala and Kleban [8] in particular. Critical for this 

^ classification are the number theoretic and geometric properties of the triangle map, 

a type of multi-dimensional continued fraction. 
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1 Introduction 

In a new classification scheme for real numbers was given, inspired by 
ideas from statistical mechanics. The key technical tools involved continued 
fractions. In that paper, it was asked if similar classification schemes occur 
for various multi-dimensional continued fraction algorithms. The goal of this 
paper is to show that such a classification scheme does happen for pairs of real 
numbers using the triangle map [26] [1] [10] [35], a type of multi-dimensional 
continued fraction. 
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Similar to we use the thermodynamic formalism (developed by Ruelle 
|32j |33j , Sinai [36] and others in an attempt to put statistical mechanics on 
a firm mathematical foundation) to define a partition function for each pair of 
real numbers. Using the real-world interpretation of the partition function, we 
classify different pairs of reals via their critical point phenomena. 

The outline of the method is as follows. We first define the relevant partition 
function. Key will be that each pair (a, j3) of real numbers gives rise to a 
different partition function, and hence can be thought of as giving rise to a 
distinct thermodynamic system. Next we use the critical point phenomena for 
each of these pairs (a, (3) to give the needed definitions for our new classification 
scheme. In order for this to be a reasonable classification scheme, we next state 
the theorems that show different pairs (a, /3) can be distinguished. In section 
four, we put our partition functions into the rhetoric of triangle sequences, 
which is the multidimensional continued fraction that plays the role of continued 
fractions from The rest of the paper contains the number theoretic proofs 
of our main theorems. 

There are two interrelated justifications for this paper. First, it shows that 
analogs of the work in llj apply for at least one multidimensional continued 
fraction algorithm. This suggests that each multidimensional continued frac- 
tion has a corresponding statistical mechanical interpretation. (For more on the 
general theory of multidimensional continued fractions, see Schweiger's Multi- 
dimensional Continued Fractions |34j ; also, many of these different algorithms 
have been put into a common framework j5] , using the triangle map as a start- 
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ing point.) Second, the proofs of these generalizations will not be obvious. The 
proofs of the two main theorems for this paper use key technical aspects of the 
triangle map, in particular its underlying geometric motivation, which is why 
the proofs of these generalizations are not just exercises. 

This current paper and its predecessor [11) are certainly not the first to 
apply statistical mechanical ideas to number theoretic questions. For example, 
there is also the transfer operator method, applied primarily to the Gauss map, 
which allows, in a natural way, tools from functional analysis to be used. This 
was pioneered by Mayer (see his |23| for a survey), and nontrivially extended 
by Prellberg [27], by Prellberg and Slawny [29], by Isola [14] and recently by 
Esposti, Isola and Knauf [B]. A good introduction to this work is in chapter 
nine of Hensley [T3] . 

Other links of statistical mechanics to number theory include the work of 
Knauf [TO] [2U] EP, of Mendes France and Tenenbaum [M] [55], of Guerra and 
Knauf [12] , of Contucci and Knauf [4 , of Fiala, Kleban and Ozliik [9] , of Kleban 
and Ozliik [16 , of Prellberg, Fiala and Kleban [25], of Feigenbaum, Procaccia 
and Tel [7], of Kallies, Ozliik, Peter and Snyder [TSJ and others. 

2 The partition function 
2.1 The General Set-up 

Motivated by statistical mechanics, we need to find an appropriate partition 
function. In section 3.1 in [11] , a general framework for developing number 
theoretic partition functions was given. Here we will simply state what our 
desired partition function should be. 
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For matrices A = (ay) and B = (bij), the Hilbert-Schmidt product (which 
is also called the Hadamard product) is 

A* B = Tr(AB T ) = ^ a H b H- 

l<i,j<n 

For example, thinking of a 3 x 3 matrix as an element of K 9 , then A * B = 
Tr(AB T ) is simply the dot product of the two vectors. 
For any (a, (3) £ R 2 , define 
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and define 



In this paper, our state space will be, for each positive integer N, 

Sat = {(<ti, ■ ■ ■ ,ctjv) : <Ji = or 1}. 
For each iV-tuple / = (cti, . . . , a N ) £ § N , set 



A — A ai A a2 ■ ■ ■ A aN . 



Definition 2.1. The partition function Z^ (a, /3, s) is 

Z N (a,[3,s)= J2 * 



By direct calculation, we have 

Z N (a,f3,s)= ]T 
ies N 
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If wc denote the third column of the matrix A 1 as 





then 



Z N (a,j3,s) = ^2 



, s lX {I) + ay{I) + pz{I)\ s - 
Thus the terms in Zn that matter most are those A 1 whose third column is 
close to being perpendicular to the vector (1 a /3 ). The behavior of Zn 
is a measure of how good of a Diophantine approximation we can achieve by 
iterations of the matrices Aq and A\ . This will be made clear in section [4] in 
our discussion of the triangle sequence. 

As discussed in if we follow the motivation and inspiration from ther- 
modynamics, we expect the partition function to pack a tremendous amount of 
information. Also, on a much more minor point, the variable s should represent 
the inverse of temperature (though in this number theoretic system, there is of 
course no true notion of temperature for numbers. 
The free energy is the function 

= lim iog(g^M 

JK ' JV-S-oo N ' 

when the limit exists. From physics, the values of s where f(s) is non-analytic 
should be important. In real world examples, it is believed that it is at these 
points where critical phenomena should occur. We will distinguish different 
pairs of real numbers based on how their corresponding number-theoretic free 
energies behave. 
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3 Classifying pairs of reals via free energy 

With the notation above, we have for each pair of real numbers (a,f3) and 
each positive integer TV the partition function Z N (a, (3, s). 

Definition 3.1. A pair of real numbers {a, (3) has a k-free energy limit if there 
is a number s c such that 



exists for all s > s c 

The subscript 'c' is used to suggest critical point. For k = 1, this is a number- 
theoretic version of the free energy of the system. By an abuse of notation we 
will also say that (a,f3) has a /(iV)-free energy limit, for an increasing function 



exists for all all s > s c . 

This will only be a reasonable classification scheme if we can show the fol- 
lowing two theorems. 

Theorem 3.2. For any positive real number k, there exists a pair of real num- 
bers (a, ft) that does not have a k-free energy limit, for any value of s. 

Theorem 3.3. Let (a,/3) be a pair of real numbers such that there is a positive 
constant C and constant d > 2 with 



lim 

N->oo 



\og(Z N (a,(3,s)) 
sN k 



f(N) if 



lim 

N-too 



\og(Z N (a,(3,s)) 
sf(N) 



1 



< \p + aq + (3r\ 
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for all relatively prime integers p, q and r, with b the maximum of \p\, \q\ and 
\r\.. Then {a,j3) has a k-free energy limit for s > 2, for any k > 1, and in fact, 
the k-free energy limit is zero. 

From Theorem 6.4 in [3J, the above yields 

Corollary 3.4. All pairs (a, /3) of algebraic numbers, each of which has degree 
at least three and for which l,a,/3 are linearly dependent over Q, have k-free 
energy limits equal to zero, for any k > 1. 

In particular, 

Corollary 3.5. For any algebraic number a of degree at least three, the pair 
(a, a 2 ) has k-free energy limits equal to zero, for any k > 1. 

4 The triangle sequence 

4.1 Review of the triangle sequence 

This is a rapid fire overview of the triangle sequence [TU] pQ [3S] . 
Partition the triangle A = {(x, y) : 1 > x > y > 0} into disjoint triangles 

A fe = {{x, y) 6 A : 1 - x - ky > > 1 - x - (fe + l)y}, 
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For a pair of real numbers (a, j3) € define 



T(a,f3) 



a 



a 



The triangle sequence for a pair (a, /3) is the infinite sequence of nonnegative 
integers (an, ai, 02, . . .), where T k (a, f3) G A afc . The map T is a one-to-one onto 
map from each subtriangle to the original A. In [55], Messaoudi, Nogueira 
and Schweiger showed that this is an ergodic map. The triangle sequence is said 
to terminate at step k if T h (a, /3) lands on the interval {(t,0) : < t < 1}. 

Now to start linking the triangle sequence with our partition function. An- 
other way of thinking about triangle sequences is as a method for producing 
integer lattice vectors in space that approximate the plane x + ay + j3z = 0. Set 

C_ 3 = 

If the triangle sequence for (a, /3) is (an, 01, 02, . . .), set 




Cfe — Cfs-3 — Ck-2 — dkCk-l- 
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The triangle sequence can be defined in terms of the dot products dk = (1, a, f3) ■ 
Cfe. Assuming we know the number an, . . . , a&, then afc+i is the nonnegative 
integer such that 

dfc-2 - rffe-i - a k+ id k > > <4_ 2 - dk-i - (a/c+i + 

Then 

For each sequence (an, Oi, 02, . . . , a*,), there is a subtriangle A(oo, • • • , Ofc) of 
A such that TC° +1 ) : A (an, . . . , a^) — > A is one-to-one and onto. In [T], the 
vertices of this triangle were determined. For this, we will need a little notation. 
Recalling the notation from [I], for vectors 

T = 

with a, d, a + d 0, define 

f 



b+e 
a+d 
c+f 
a+d 




and further, define 

T+S = T+~S = 

(Such a sum is called a Farey sum.) Further, continuing with this type of 
notation, we define 

a \ I d \ ( g \ ( b+e+h \ 

C \ f \ f V -+d+9 J 



As shown in pQ, we can now cleanly describe the vertices for the triangle 
A(ao, • • ■ ,a n ). Set 

X n = I Dn = C n X C n +i. 
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Then the vertices for the triangle A(an, . . . , a n ) are X n —i, X n and X n +X n -2- 
At a critical stage for our argument, we will be using that 

A(a,X k -i,X k ) c A(X k +ak+iX k - 1 +X k - 2 ,X k - 1 ,X k ), 

which will mean that area A (a, X k -i,X k ) < area A(X k +a k +iX k -i+X k -2, -Xfe-ij X k ), 
as seen in 




4.2 The triangle sequence and the A and A\ 

Here we provide the key link between the geometry of the triangle sequence 
and the terms for our partition function. Let (a, /3) € A have triangle sequence 
(ai, <i2, . . .). For each k, set 

7V fc = (ai + 1) + . . . (o* + 1). 
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Then, as seen in [T], the matrix A^ 1 A$A a2 Ao ■ ■ ■ A ak Ao will have columns CN k -2, Cm 
Thus for the sequence I = (l ai , 0, 1 Q2 , 0, . . . , l a " , 0), we have 

M* A 1 = d Nk . 

5 Proof of Theorem 13.21 

Using the notation of the previous two sections, we start with 
Lemma 5.1. For all k, 

1(1 a (3))-C k \< . 

\x k +i\ 

Proof. We know that 

a e A(X k +X k _ 2 , X k -i, X k ). 
The definition of X k+ \ and of a k +i is that 

but that 

( a ) £ A(X k + {a k+1 + l)X k -!+X k - 2 , X k - U X k ). 
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|(1 a P)-C k )\ = |(1 a )-(X k _ 1 xX k )\ 

= |det(( 1 a p ),X k _ 1 ,X k )\ 

= |x fc _ix fe areaA(( a P ) , Xk-i, X k )\ 

< \x k -iX k axe&A(X k +a k+ iX k -i+X k - 2 , X k -i,X k )\ 
detiXk + Xk^Xk^Xk) 



Xk~ix k {xk + a k +ix k -i + x k _ 2 ) 
1 



\x k + a k+ ix k -i + x fe _ 2 | 
1 



\x k +i\' 

as desired. 

Here we are using the following two facts. First, if 
A = 

the area of the triangle with vertices 
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ai3 
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Second, we have 



det(A) 

a llOl2 a 13 



det(X k + X k - 2 ,X k -i,X k ) = 1, 

which is stems from the fact that the matrix (X k + X k _2i X k ^i, X k ) is the 
product of matrices of determinant one. 



□ 
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We are now ready to prove Theorem 3.2 namely given any positive real 
number k, there exists a pair of real numbers (a, (3) that does not have a fc-free 
energy limit, for any value of s. 

Let f(n) be a function that strictly increases to infinity. Choose a sequence 
of integers 0,1,0,3,... inductively by setting 

„ , , > p f(m+l)(a 1 +...a m +m) k 

with the initial a% any positive integer greater than two. 

Choose a pair (a, 0) so that the pair's triangle sequence is (0%, et 2 , . . .). Recall 
that we set N m = a± + 02 H + a m + m. 

Consider Z Nm (a, /?, s). We have the sequence I = (l ai , 0, 1 Q2 , 0, . . . , l a ™ , 0) 
being part of what is summed over for the partition function Zpf m (a, ft, s) . Thus, 
for this /, we have 

1 _ 1 

|(1 a J) )-C N J" ~ \M*A Z \' ^ \M * A'\ 



|x m+ i| s < Z Nm (a,0,s). 



\og\x m+1 \ s 

lim r = OO 

m— >oo ]S[ 



Thus we have 

We will show that 

for any value of s > 0. 
We have 

Thus if m is even we have 

%m+l > O'm+lO'm-l ' ' ' a l 
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In particular, we have 



X m +1 > <hn+l a m-l " " * a 2- 



Xm+1 > flm+1 



Then 



log|a; m+ i| s log |a m+ i| s 

urn , > hm 



m->oo (ai H a m + m) k 

S log \ e f(m+l)(ai+...a m +m) k | 



> lim 



(ai + ■ ■ ■ a m + m) k 
sf(m + l)(ai + . . . a m + ?7i-) fe 



= lim 

tn-s-oo (ai 4- • • • a m 4- m) K 

= oo, 

finishing the proof of Theorem |3.2| 

6 Proof of Theorem 13.31 

By assumption, (a, j3) is a pair of real numbers such that there is a positive 
constant C and constant d > 2 with 

-^j < |p4-ag4-/3r| 

for all relatively prime integers p, q and r, with 6 the maximum of |g| and 
|r|.. We must show that (a, 0) has a k-free energy limit for s > 2, for any fc > 1. 
In fact, we will see that the k-free energy limit is zero. 

By inverting and raising both sides of the above inequality by s, we have 

1 



\p + aq + (3r\ : 



< C s b 



s isd 
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In this case 



Z N (a,f3,s) = ^2 



ies t 



\x{I) + ay{I) + pz(I)\ i 



< C s (max(|x(I),y(I),z(I)) s 



Hence we need to control the growth rates of the \y{I)\ and \z(I)\ in terms 

of the integer N. We will find fairly crude but workable bounds. First off change 
the minus ones in the matrices Aq and A\ to ones. Then the growth rates of the 
|x(J)|, \y(I)\ and \z(I)\ will be bounded by the N*h Fibanacci number, which is 
well known to be 



V5 ' 



where 



Thus 



1 + V5 , 

r, ! ^ = n ■ 



sd 

Z N (a,p,s) < 2 N C s F s jf < 2 N C S ' 



V5 



Hence 

],,,, [ 2 N C S 1 — % " 

hm l0g(Zj ! (a;5)) < lim 



N-Hx (3N k ~ N->oo sN k 

N log(2) + s log(C) - sd log("v/5) + sdN log(0) 

N^oo sN k 

= 0, 



for K > 1, as desired. 
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7 Conclusion 

To some extent, this paper is the second (with [11] the first) in a possible 
series of paper using the thermodynamic formalism to understand more about 
the structure of the real numbers. There are many questions left. 

Probably the most straightforward, but still non-trivial, ones are to try to 
mimic the results in this paper for other types of multidimensional continued 
fractions. These analogs would involve choosing different matrices Aq and Ai. 

Our partition function Zw(a, ft, s) depended on the choice of the three-by- 
three matrix M. Different choices of M give rise to different thermodynamics. 
Certainly choosing different M's in the two-by-two case led to the differences 
in the work of Knauf PHQ1] HD] and of Fiala, Kleban and Ozluk [g][H] 

versus the work in [llj . Are there analogs here? 

We have shown that some pairs of real numbers have 1-free energy limits 
while others do not. This too is just a beginning. What types of limits are 
possible? Are there pairs (a, (3) for which the sequence lo s(M n 'ft s )) h as anv 
possible limit behavior? For example, can we rig (a, /3) so that any number can 
be the limit of the sequence. In fact, we should be able to find such sequences 
with two accumulation points, three accumulation points, etc. All of these 
should provide information about the pair of real numbers a and /3. 

In |11] , continued fractions were critical. Since an amazing amount is known 
about continued fractions, it is not surprising that there are more results in 
|llj . For example, in that paper it was shown that e — 1 has a one-free energy 
limit. Are there similar such results possible for concrete pairs of reals of the 
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thermodynamic systems of this paper. These types of questions are to some 
extent the most intriguing. 

Finally, there is the question of putting these results into the language of 
transfer operators. (See [3T], [35], [2] for general references.) 
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